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Utilizing a novel and exact inversion scheme, we determine the complex linear conductivity <j(u) 
from the linear magnetic ac-susceptibility which has been measured from 3 mHz to 50 MHz in fields 
between 0.4 T and 4T applied parallel to the c-axis of a 250 nm thin disk. The frequency derivative 
of the phase a" /a' and the dynamical scaling of a(u) above and below T g (B) provide clear evidence 
for a continuous phase transition at T g to a generic superconducting state. Based on the vortex- 
glass scaling model, the resulting critical exponents v and z are close to those frequently obtained 
on films by other means and associated with an 'isotropic' vortex glass. The field effect on a(w) can 
be related to the increase of the glass coherence length, £ 9 ~ B. 



PACS numbers: 74.25 Nf, 74.60 Ge, 74.72 Bk, 74.76 Bz 



I. INTRODUCTION 

The widely observed scaling feature of the electrical 
conductivity in the presence of a transverse magnetic field 
B and the universality of the critical exponents obtaij 
in such scaling analyses on thin YBa2Cu307-filmsB~ 
serve as the best evidences for the existence of a continu- 
ous thermodynamic phase transition to a genuine super- 
conducting phase occuring at some transition line T g {B). 
Considering the interplay between the elastic properties 
of the vortex lattice and the disordering potential ex- 
erted on the vortices by density fluctuations of point- 
like pinning centers on a microscopic level, the possibil- 
ity of a glassy vortex state with infinite barriers for vor- 
tex motion, U(j — > 0) and, hence, zero linear resistance, 
p(J ~ v 0) = 0, was argued previously for 3-dimensional 
lattices (see e.g. Ref. 12 and references therein). How- 
ever, a theoretical nronf. for the hypothesized second or- 
der phase transitionEjO to this vortex glass (VG) which 
is characterized by both a singular correlation length and 
diverging relaxation time of the order parameter fluctu- 
ations, 



S g (T,B)=£ g (B)\l-T/T g 



T 9 (T,B)^r g (B)\l-T/T g 



(la) 



(lb) 



is still lacking. Hence, also the analogy between the tran- 
sition to the Meissner-phase at B = (relaxational 3d- 
XY-model) and to the VG-phase at B ^ defining the 
scaling property of the conductivity does not have a firm 
basis yet. Within this model, magnitude and phase angle 
are predicted to obey the following formai-J 



\a(u > )\ = fS ± (ur).. 



arctan 



(2a) 
(2b) 



where S±(x) and P±(x) represent homogeneous scaling 
functions above and below T g , which describe the change 
from pure ohmic behaviour (a" = 0) for T T g to pure 
screening (a' = 0) at T <C T g . While their limiting be- 
haviours are known, the detailed shapes of these scaling 
functions have not yet been worked out. 

This somewhat unsatisfactory theoretical situation 
contrasts to the ample evidence for scaling and universal- 
ity obtained on thin YBCO-films. Most of these analy- 
ses were based on the nonlinear dc-conductivity_CT(j)Ej1j, 
a recent onp-.an the non-linear Hall resistances, and in 
two studiestTJO the linear dynamic conductivity at T g 
has been investigated. Within the given error margins, 
almost all of the published critical exponents can be sum- 
marized by the r.m.s. averages v = 1.7(1) and z = 5.5(5). 
This represents a nice example for universality, since 
most films have been prepared under different conditions 
and do deviate from each other in microscopic composi- 
tion and structure, normal resistivity, transition temper- 
ature and width. Moreover, the exponents appear to be 
insensitive against variation of the vortex density B/$>o 
and direction of the field with respect to the c-axiso. 

On the other hand, notable exceptions fr,om this uni- 
versality have been seen. Dekker et al.0 observed a 
crossover to two-dimensional-critical behaviour when the 
film thickness is reduced to about 16 A, i.e. to the 
thermal correlation length I-V data from micro- 

bridges (80 x 10/im 2 ) patterned from epitaxial YBCO- 
films have been interpretedt£l in terms of thermal acti- 
vation across Josephson-barriers using the Ambegaokar- 
Halperin modeLj. For YBCO-crystals distinctly differ- 
ent exponents have been reported (see e.g. Ref. 17) indi- 
cating the possibility of another universality class for the 
transition to true superconductivity due to the presence 
of a twin boundary network. 

This paper reports a first detailed experimental study 
of the scaling behaviour of the linear dynamic conduc- 
tivity ct(lo) of a thin YBCO-film in a wide frequency 
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range covering more than ten decades from 3 mHz to 
50 MHz. This compares to the existing measpseaients of 
a(co) on films between 0.5 MHz and 500MHzNo which, 
moreover, focused only on the critical behaviour of <r(u) 
at T = T g . Using experimentally well defined scaling 
functions for <r(w), we will also examine the not yet ex- 
plained effect of the magnetic field on the vortex-glass 
scaling, at least on an empirical level. Recently, within 
the mean field approximation!!!!, the scaling behaviour 
of <t(lu) has been studied on the basis of the Edwards- 
Anderson model. 

Our experimental method employs a contact-free ac- 
cess to the dynamic conductivity of a thin film. We mea- 
sure the linear dynamic susceptibility x(cj) parallel to 
both the applied dc-field and the c-axis, which is normal 
to the film plane (see section II). Using an exact inver- 
sion schemeEj outlined in section III, we then evaluate 
from x( w ) the complex linear conductivity a(uj) in the 
film plane. Section IV presents the scaling analysis of 
er(cj) and determines the relevant results, i.e. the criti- 
cal exponents, the glass line, and the field effect on the 
scaling functions. The results are discussed in section V. 



II. EXPERIMENTAL 

The 250 nm thick disk (R=1.5mm) under investiga- 
tion has been cut from a film, epitaxially grown parallel 
to the (OOl)-axis of a SrTi03-substrate by pulsed laser de- 
position being described in detail elsewhereo. AFM im- 
ages show a rather flat surface from which some isolated 
hills of 1 /Ltm average diameter and up to 0.2 /zm height 
grew up. The dc-resistance at 100 K is about 100/if2cm 
and the transition temperature is T c = 91.4 K with a 
10:90% width of less than 1 K. From the non-linear ac- 
susceptibility critical current densities of 2 • 10 6 A/cm 2 
and 10 7 A/cm 2 were obtained at 77 K and 20 K , respec- 
tively. 

In order to measure the linear susceptibility, the sam- 
ple, a thin circular disk (R = 1.5 mm) was inserted in 
well balanced mutual inductances, the excitation am- 
plitudes of which were always controlled to obtain the 
linear response. We use two different detection tech- 
niques: between 3 mHz and 1kHz a commercial SQUID 
ac-susceptometer (Quantum Design MPMS2) and from 
10 Hz to 50 MHz conventional lock-in amplifiers (PAR 
5301 and PAR 5202). At fixed external fields, provided 
by superconducting magnets in the persistent mode, two 
measuring schemes have been used. Most of the data 
were recorded by a quasi static variation of the tempera- 
ture at constant frequencies so that the same signals were 
obtained during heating and warming of the sample. Sev- 
eral measuring curves are reproduced by Fig. 1. At some 
interesting points in the B-T phase diagram the tempera- 
ture was stabilized to investigate the frequency variation 
towards the low frequency end of our spectrum, where 
the measuring time for a magnetic moment of 10 -7 emu 



to a nominal accuracy of < 10 -8 emu increased up to 12 
hours. In the entire frequency and temperature range of 
interest here, the complex output of the mutual induc- 
tances has been carefully calibrated against the known 
signals of para- and ferromagnetic samples (Pd, EuO, 
GdCi3- 7H2O) of the appropriate volumes. 



III. AC-CONDUCTIVITY 

To determine the dynamic susceptibility x( w ) = 
(l/2Vh ac ) J* x j(r, to) d 3 r of a thin disk of thickness d 
in a perpendicular ac field h ac e luJt one has to solve an 
integral equation which describes nonlocal diffusion of 
the normalized in-plane sheet current density J(r, z) = 
j(r, z)d/h ac . In the following, complete penetration \5\ > 
d is assumed, equivalent to lu < 2/ /j, d 2 \a ac \ with the skin 
depth 5 = (2/ hqujct) 1 / 2 . This integral equation follows by 
inserting the nonlocal relation between J(r, z) and the 
perpendicular field h x (r,z) it generates (Ampere's law) 
into the induction law V x e = — b using the appropriate 
material laws j = a ac e and b = /ioh. For a disk of radius 
R (r < R) this integral equation reads 



J(r, lo) — w nr + / Q(r, u) J(u, lo) du 
Jo 



where 



iuRdnQ 
w = — cr ac (u) 

Z7T 



(3) 



(3a) 



In the Ohmic case, a ac = <r is real, thus w is purely imag- 
inary, whereas in the Meissner state or for rigid pinning, 
a ac is purely imaginary and thus w is, real and frequency- 
independent. The integral kernel ia!3 



Q(r, u) 



r/u 



E(k) K{k) 



1 — v 1 + v 



v dv . 



(3b) 



Here E(k) and K(k) are complete elliptic integrals and 
k 2 = 4?;/(l + v) 2 . Expanding J(r,uj) in terms of the 
eigenfunctions f n {v) (Fig. 2) of the eigenvalue problem 



fn(r) = -A„ / Q(r, u) f n (u) du , 
Jo 



(4) 



one gets for the planar current density J(r, lu) — 
Ln a «( w )/™( r ) witn 



a n (u>) 



ir wb„ 



1 + w/A n 
yielding the susceptibility 



f n {r)r 2 dr. 



X 



The orthonormality, f Q fn(r)f m (r)r dr = S mn of the 
eigenfunctions was used. The eigenvalues are A„ rj 
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n — 0.123 for all n = 1,2,3,.... Finally, one obtains 
in the form 

x (u))=tt[ J(r,uj)r 2 dr = -wV c ™l^~ n ) ( 6 ) 

JO n A " + w 

with c„ = (3-7r 2 /8)6^A^. This quite general result ex- 
presses as an infinite sum (n = 1... oo) of terms 
which have first order poles at real positions A„ in the 
complex w plane, or at complex positions in the lj plane. 
In the Ohmic case, the poles all lie on the positive imagi- 
nary axis of the uj plane. If a ac (uj) = |<r ac |e lP is complex, 
the poles have the phase angles tt/2 — P(u->). 

Next we show how the amplitudes c„ and eigenvalues 
A„ may be calculated and how the infinite sum (||) can 
be approximated with high precision by a finite sum of 
N terms (n — 1 ... N). After discretizing the continuous 
variables r and u as described in Ref. 19, the eigenvalue 
problem (|4|) is equivalent to the diagonalization of an 
N x N matrix Qij defined by 

n nf ^ n , 0.92363^- 

Qi& = Wj Q [Ui, Uj) , Qjj = Wj In — — — — . 

Here the Ui — u[Xi) span a nonequidistant grid ob- 
tained by inserting equidistant x% = (i — i) /N (i — 
1...N) into an appropriate substitution function, e.g. 
u(x) = (75x — 50a; 3 + 15x 5 )/40 which yields a weight 
function w(x) — u'(x) = (15/8)(1 — x 2 ) 2 and weights 
Wi = w(xi)/N. This grid, and the diagonal terms 
Qjj, are chosen to maximize the accuracy of the sums, 
which approximate integrals like J Q L /(r)dr rs J^/jitf, or 

JoQ( r i> u )f( u )du ~ J2Qijfj> see R- ei - 19. The high ac- 
curacy of this grid (error ~ N~ 3 ) is achieved by the van- 
ishing of the weight function w(x) at the edges of the 
disk, x = u = r = 1, where the arbitary function f(r) 
may be peaked or may have a logarithmically diverging 
slope. In fact, the eigenfunctions f n (r) and the integral 
kernel Q(r,u) exhibit such weak singular behaviour at 
the edge r = 1, which is hardly seen on the linear scale 
of Fig. 2. 

The matrix Qij has exactly N eigenvalues A n and 
N eigenfunctions /„. Inserting these into (||) one gets 
the coefficients b n and the A" amplitudes c n entering 
the sum (g) for x( w ); see Table I. The finite sums ob- 
tained in this way present the exact solution of the dis- 
cretized problem and are very good approximations Jx> 
the original problem, the magnetic response of the disklij. 
The sum rules tf n = 1/4, £„^A„ = 8/3tt 2 , and 
Sn^n^n 1 = 4/15 (see Ref. 19) are satisfied with high 
precision even for small A^ rj 10. However, while for 
A^ — > oo the b 2 n « 8/7T 4 n 3 and c ra w 37r 2 jn decrease mono- 
tonically with n and the A„ w n — 0.123 are equidistant, 
finite N < oo yields nonmonotonic amplitudes c„ and 
nonequidistant eigenvalues A„ (Table 1) in order to com- 
ply with these sum rules and fit the exact infinite sum 
(6) for x(w) in a large range of to. 



In order to invert the measured ac-susceptibilities 
x'(k>) — ix"(w) into ac-conductivities <j'(oj) + ia"(u>) by 
means of the analytical expression (^), we used a numer- 
ical inversion (searching) procedure. The choices N = 30 
(Table I) and N = 50 yielded identical results within the 
accuracy of our data. 

As an example, Fig. 3 displays the phase angle and 
the modulus of cr(ui) evaluated from x( w ) at 0.4 T (see 
Fig. 1(a)). Even slightly below T c these data exhibit a 
frequency dependence which cannot be explained within 
the proposed extended models of thermally activated 
flux flow (TAFF)cll. The most striking difference oc- 
curs in the phase P, which according to TAFF contin- 
uously increases from zero to n/2 when passing uj rj 
T _1 exp(— U/T) and then decreases in the same manner 
to zero above the much higher frequency u ~ t^ 1 due 
to losses originating from the vortex motion within their 
pinning wells. This is in stark contrast to the present be- 
haviour (Fig. 3) where, with increasing temperature, the 
increase of the phase P (dP/ doj > 0) comes to a halt be- 
low P = tt/2 and is followed by a region with dP/dui < 
towards P — tt/2. Hence our data indicate the onset of 
collective effects in the fluid vortex phase near T c , which 
are not included in the TAFF model. 



IV. SCALING ANALYSES 

When comparing Figs. 1 and 3, the interesting feature 
emerges that xi 10 ) does not display any direct signal of 
a thermodynamic phase transition, whereas the phase 
arctan(cr"/c') of the ac-conductivity provides strong ev- 
idence for such a transition: its frequency derivative 
changes sign at some finite value P below tt/2 demanded 
by the-continuity of the scaling function P±(x) (Eq. (2b)) 
at Tgtil. The physical reason behind this characteris- 
tic feature can be associated with the coexistence of in- 
finitely large clusters (£ s — > oo) of superconducting glassy 
and a metallic (liquid) vortex phase at T g . With increas- 
ing temperature, the life time and size of the glassy fluc- 
tuations decrease so that within larger observation times 
l/u> they cannot contribute to the screening component 
a", whereas towards lower temperatures the reverse is 
true: within increasing l/u> the lossy contributions a 1 
arising from vortex loop excitations die out at the ex- 
pense of the stable glassy phase. 

Basing on this signature for the continuous phase tran- 
sition at T g , the next step is to check the scaling property 
predicted for the phase angle (Eq. (2b)) in terms of the 
variable lot ~ u\l — T/T g \~ vz by a proper adjustment of 
the exponent vz. In fact, Fig. 4(a) shows that by taking 
vz = 9.35 the data scale very nicely over more than 14 
decades in the scaled frequency lot. Together with the 
value z = 5.5(5) following from the phase angle at the 
transition P(T g ) = (1 — 1/z)tt/2^3 also the exponent v 
is determined, v = 1.7(1). The 4T-data also shown in 
Fig. 4(a) prove that both exponents are independent of 
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the magnetic field, at least between 0.4 T and 4T stud- 
ied here. Moreover, these results are identical with the 
averages for v and z, defined by the numerous scaling 
analyses of the nonlinear conductivity of thin YBCO- 
films as outlined in the introduction. This constitutes an 
additional independent confirmation of the universality 
feature for YBCO films. Note that this transition from 
'metallic' flux-flow behaviour (P = 0) to true supercon- 
ductivity (P = 7r/2) occurs in rather narrow temperature 
intervalls. For 0.4 T and 4T, we find the <r(a;)-data from 
the intervalls ±0.9 K and ±3.0 K around T g to collapse 
on the scaling functions. These intervalls correspond to 
reduced temperatures |1 — T/T g \ smaller than 10~ 2 and 
3 • 10~ 2 as expected for a critical phenomenon. Inter- 
estingly, they are in almost quantitave .agreement with 
the critical regions quoted by Reed et alE3 on a twinned 
crystal. 

Additional support for the scaling property of u(w) 
arises from the modulus of the conductivity shown in 
Fig. 4(b) for the same two vortex densities. Here the 
natural scaling variable is the dc-conductivity above T g , 
°s(0) ~ T s/£s ~ ll-T/Tg]-^- 1 ^. Due to the validity of 
the Kramers-Kronig relations for cr(u>), no additional pa- 
rameter is required to determine the scaling according to 
Eq. (2a). Using essentially the same glass temperatures 
and critical exponents, excellent scaling of all data from 
the indicated temperature intervals above and below T g 
is found. Obviously, they collapse on two branches defin- 
ing the scaling functions S±(x). Again, their limiting 
cases S+(0) and S_(0) represent the pure lossy (|ct| = a') 
and pure screening (|<r| ~ a" ~ w _1 ) responses, while 
at T g one finds the algebraic behaviour \a\ ~ oj~ 1+1 I z 
predicted for the dynamical criticality. Principally, these 
results for \a\ can be also considered as a successful check 
of the measuring and new inversion procedures presented 
here. 

The inset to Fig. 5(b) shows the phase diagram defined 
by these scaling analyses. The transition temperatures 
define the glass line separating the normal conducting 
vortex fluid from the superconducting vortex glass, which 
can be described by the power law: 

Bg (T)=B g (O)(l-0 , (7) 

where f3 = 1.5(2) and B g (0) = 85(5)T. 

Looking at the scaling functions S± and P± exempli- 
fied in Fig. 4 for 0.4 T and 4 T, one notes that their shapes 
appear rather similar to each other. This strongly sug- 
gests to perform an additional scaling procedure in a phe- 
nomenological manner. We simply assign suitable pow- 
ers of the applied field B to the scaled frequency and 
conductivity. In fact, Fig. 5 demonstrates nicely that 
for both, phase and modulus of cr(u>), all experimental 
scaling functions fall on single curves over 15 decades of 
the scaled frequency if we assume that the lifetime r g of 
the glassy fluctutions increases with 5 5 - 5 ( 10 ) and the dc- 
conductivity, <t(0) ~ T g /^ g (Eq. (2a)), grows proportional 



to £} 4 - 5 ( - 5 ). Interestingly, both results can be related to 
the same physical origin, namely to the field dependence 
of the vortex-glass correlation length. We obtain the field 
dependence of the principal scaling variable in the follow- 
ing form 

UT,B)=UBo) ^ / / T ^°}\^ , (8) 

where 7=0.59(8) and Bq denotes an arbitrary normal- 
ization field. This finding implies that the shape of the 
scaling functions does not depend explicitly on the mag- 
netic field, but that the field effect on the scaling can be 
fully accounted for by the B dependence of the glass cor- 
relation length £ g (T,B). Moreover, if we follow Ref. 17 
to express the field dependence of £ s in terms of the dis- 
tance between T g and T C 2, i.e. \T C 2 — T g \ ~ T g B@ , with 
/3' = 1.6(2) for the present glass line, we find within the 

given error margins, £ /v (T,B) ~ (T c2 - T g )/\T=T g \. 
This is fully consistent with the previous resultliU ob- 
tained on a single crystal. 

V. DISCUSSION 

The main information obtained from scaling analy- 
ses are the critical exponents, which are determined by 
the rather general features of the phase transition, i.e. 
the spatial dimensionality, the symmetry of the order- 
parameter and the type of disorder. The present values 
for the exponents of the correlation length and the relax- 
ation time, v — 1-7(1) and z — 5.5(5), respectively, agree 
nicely with Hiftst of the existing values determined for 
YBCO-filmsEJ eI from non-linear I-V characteristics. On 
the experimental side, this provides a strong support for 
the ac-susceptibility technique and the evaluation proce- 
dure applied here for the first time on a superconducting 
thin film. For bulk YBCO-crystals such an exact trans- 
formation of the complex susceptibility to the dynamic 
conductivity using conventional schemes has been per- 
formed recentlyO, while another recent scaling analysis 
was performed for the modulus.-af resting on the 

assumption |x(w)| — a/|A ac (w)|E3. Of course, this as- 
sumption is only valid for A oc being large compared to 
the sample dimension and should be limited to small sus- 
ceptibilities, i.e. to just the onset of screening. 

The present values of the exponents z'.-and z are signifi- 
cantly larger than the mean-field valuesc£l for the vortex- 
glass (VG) v — 0.5 and z = 4, indicating the pres- 
ence of strong critical fluctuations. For the case of a 
3-dimensional isotropic VG one expects from analogy 
with spin-glassesEfl larger numbers, e.g. z = 2(2 — rf) 
with t] < dl3, however, exact values, are not known 
yet. Monte-Carlo simulations of thp-gauge-glass model 
revealed v = 1.3(4) and z = 4.7(7)EC but the applica- 
bility to the vortex state in real disordered materials like 
YBCO-films has not been proven to date. 
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On the other hand, the results v = z = 3.1(3) from 
scaling .analyses of u(ui) for B \\ c of twinned YBCO- 
crystalsOO indicate the existence of a distinctly differ- 
ent universality class for the transition to a supercon- 
ducting vortex state. It has been suggested to associate 
this class with correlated pinning along the boundaries-of 
twin colonies (|| c), where single twin planes terminateE3. 
In this case of anisotropic pinning, one has to distin- 
guish between two correlation lengths £|| and £j_ of glassy 
fluctuations parallel and perpendicular to the preferred 
c-axis. Then the measured exponents can be explained 
by the assumption £|| ~ Q_. These anisotropic fluctu- 
ations are considered by the so-called Bpse-glass model 
for a vortex-state with columnar pinsE3, for which by 
analogy with the dirty boson problem exists now good 
evidence for a phase-transition to genuine superconduc- 
tivity at a finite temperature Tb g (B)c3. Similar evidence 
is still lacking for the isotropic VG-model, where the vor- 
tices are pinned by fluctuations of pointlike defects of 
extension r$ <C £t- Interestingly for the YBCO-crystal, 
the isotropic values of_the exponents v and z have been 
recovered for B _!_ c E3, i.e. for fields perpendicular to 
the columnar pinning centers, which means that these 
line defects are ineffective in this geometry. In contrast, 
when passing from B _L c to B || c for thin YBCO-films, 
no change of the isotropic exponents occursD, which re- 
veals an absence of columnar pinning. For our film such 
a conclusion for isotropic pinning is supported by the 
fact that twin colonies could not be detected down to a 
length-scale of 0.8 ^m. 

Another universal property of VG seems to be the ef- 
fect of the magnetic field on the transition. The field de- 
pendence of the glass-line could frequently be explained 
by the power law, Eq. (7), with exponents similar to the 
one obtained here, f3 — 1.5(2). Usually, f3 — 2i/q is re- 
lated to the exponent of the thermal correlation length 
of the homogeneous statdH £ T ~ (1 - T/T c )- U °. On 
the other hand, a detailed understanding of the non- 
universal amplitude-|£> g (0) , i.e. its dependence on the 
pinning mechanismli-3 has not yet been achieved. 

Much less evidence exists for the field dependence of 
the dynamic scaling functions S± and P± themselves. 
Our analysis between 0.4T and 4T, shown in Fig. 5, sug- 
gests that their shape is not influenced by the vortex 
density, but that the entire effect on the dynamic con- 
ductivity near T g can be accounted for by that of the 
VG correlation length. Our result, Eq. (8), implies, that 
the amplitudes in (Eq. (1)), £, g (B) and r g (B), increase 
approximately linearly with B (due to 71/ « 1.0(2)) and 
with B z , respectively. Because this finding agrees with 
that of our previous study on an YBCjQ. single crystal 
for both orientations of the applied fieldliZl, it appears to 
be another rather general property of the VG-transition. 
The most obvious consequence of £ ff ~ B is that below 
T g the squared screening length A 2 defined at low fre- 
quencies increases linearly with B: 



0,T<T g ,B) 



Mo^A 2 . 



X 2 S (T,B) 



B 



(1 - T/T g y B 



(9) 



(9a) 



with A S (1T) = 0.7/mi. 

Let us first note that this behaviour agrees with that 
of the elastic-Dinning length, Xfj = B§ /C p , introduced 
by CampbelS for a vortex lattice exposed to a pinning 
energy density C p . In this classical model, the screening 
length A s results from the elastic response of the vortex 
lattice and diverges near T C 2(B), where the lattice melts. 
According to the present results, this melting occurs at 
T g (B) and the elastic response is determined by a field- 
independent pinning energy C p . 

When we associate the screening behaviour, Eq. (9) 
with the VG-fluctuations, as suggested by the scaling 
property (Fig. 5), then the screening depth A s is related 
to the VG correlation length by the form 



At 



(10) 



where__AT = < E > o/47r/iofcs7 1 represents the thermal 
lengtltj. Inserting the experimental X s we obtain for the 
amplitude defined by Eq. (la), £ g (B) w 2.5nm x B(T). 
For the present fields, this value is in the order of the ther- 
mal coherence length at zero temperature £t(0) which is 
a natural lower bound to £ s at T = 0. This implies that, 
in the critical region, 0.002 < |1 - T/T g \ < 0.03, where 
scaling is observed here, £ g varies from about 1/im to 
100/xm. This range is plausible, since it is located be- 
tween the vortex-lattice constant, e.g. do ~ 50nm at IT 
and the radius of our disk, R = 1.5mm, which both set 
natural boundsJjor the range of the fluctuations. 

Fisher et alo argued that £, g (B) may be associated 
with the collectpe pinning length introduced by Larkin 
and Ovchinikovo, which results from the competition 
of the elastic and the pinning energies, l p w ^(cei/cp) 2 , 
where U is the mean distance between the pinning cen- 
ters. Taking £ g ~ l p for granted, our result implies 
c e i/c p ~ i? 1 / 2 which is not ampiig the predictions for 
l p in the various pinning regimesE3. Since the theories of 
collective pinning do not considder critical fluctuations 
but rather apply to T T g this discrepancy should not 
be taken as too serious. 

One interesting consequence of our finding £ 9 ~ B is 
that due to r g ~ £f ~ B 5 5 the relaxation time of the 
fluctuations is growing dramatically with field. In partic- 
ular, this leads to an increase of the dc-limit of the fluc- 
tuation conductivity, er(0) ~ ~ B 4 ' 5 as evidenced 
by Fig. 5(b). This is in sharp contrast to the conven- 
tional flux-flow behaviour, as ~ B^ 1 , and indicates that 
even slightly below T c the dc-limit of <j{u>) in the present 
YBCO-film is dominated by cooperative effects. Most 
likely, this is also true for other films with glass lines 
T g (B) sufficiently close to T c . 
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Let us finally recall the physical differeppe between the 
vortex dynamics according to the TAFFo and relatedta 
models and the highly collective VG-model. In the for- 
mer, single vortices or bundles become increasingly de- 
pinned at elevated temperatures, and this is described 
by a single activation rate r _1 (T). At fixed frequency, 
more and more vortices depin within a time w _1 and 
generate electric fields, i.e. ac-losses; this implies that 
the phase angle P = arctan(cr"/cr'), decreases gradually 
from 7r/2 to zero with increasing lot. In the VG picture, 
on the other hand, a flow of vortex bundles of extension 
£ g occurs only above the phase transition T g (B). At T g 
bundles of infinite linear extension appear, which coexist 
with other, infinite clusters of the vortex-liquid. Due to 
this coexistence of both glassy and liquid clusters with ex- 
tremely long relaxation rates, T g ~ £?, the phase angle P 
assumes a nontrivial value (1 — l/z)n/2, which lies below 
the superconducting limit P = tt/2. This limit is reached 
for lower temperatures, where the size £, g (T < T g ) of 
the vortex-liquid clusters, i.e. of the "metallic" regions 
shrinks. Since their lifetime is then strongly reduced, 
these clusters do not cause losses at low frequencies 
uo < r _1 . It is this process which gives rise to the change 
of the slope of the phase dPj dto at T g and is not contained 
in the TAFF modeM In the extended TAFF-modelN, 
a negative slope of P{u>) appears at frequencies, where 
the (viscous) damping of the pinned vortices is reached. 
In contrast to the VG-transition, the crossover from the 
low-frequency TAFF-regime (dP/du > 0) to the viscous 
damping one (dP/duj < 0) is of purely dynamic origin: 
at finite temperatures, there is no diverging time scale, 
and as long as the activation energy is large enough, the 
slope of the phase changes its sign always at P = tt/2, 
and the high-frequency limit is P = as compared to 
(1 — z~ 1 )tt/2 in the VG-model. This dynamic crossover 
to the London-screening, where \ s = Ax,, will be a matter 
of future research in the VG-state. 



VI. CONCLUSIONS 

In conclusion, we have presented an extensive study of 
the linear ac-susceptibility measured perpendicular to a 
thin circular YBCO disk, which qualitatively shows the 
transition from full penetration to full screening of the 
ac-field with increasing frequency and decreasing temper- 
ature and field. In contrast to previous investigations, 
which discuss just the shift of the maximum of x" m 
terms of the TAFF or other conventional pinning mod- 
els, we have evaluated here the full information available 
from the data for x( w > T, B) = \' — ix" ■ This evaluation 
was achieved by a novel exact inversion routine, which al- 
lows to determine the linear dynamic conductivity of the 
bulk material, u{ui). The frequency dependence of the 
phase, (t"/ct', provides clear evidence for the existence of 
a continuous phase transition to a truly superconducting 
state at a well defined temperature T g (B). Using a step- 



by-step evaluation of the dynamical scaling proposed for 
this transition it was possible to extract the critical ex- 
ponents v and z, which agree, adth those obtained on a 
large number of YBCO-filmsmlil utilizing the nonlinear 
dc-conductivity. Since these films differ in their micro- 
scopic structure as well as in their transition tempera- 
tures, the present agreement obtained with a new exper- 
imental method, strongly supports the universality which 
has been conjectured for a transition to a vortex-glasstfl 

The dynamical scaling functions for the linear conduc- 
tivity determined here for the first time in some detail, 
perhaps may serve as a challenge to the theory to work 
out a realistic model of the vortex-glass state and its 
critical behaviour. As a hint for such an ansatz one 
may consider the observed linear field dependence of the 
coherence length £ g (T, B) associated with this transi- 
tion. Other interesting questions are the crossovers from 
the critical regime to the London-screening on the low- 
temperature side and to the flux-flow limit, at T 3> T g . 
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TABLE I. Positions of the poles A n and amplitudes c n 
entering the transverse susceptibility (6) of a circular disk. 
With these numbers inserted, the finite sum (6) approximates 
the complex function \ of the complex argument w with high 
precision up to |w| « 4000. These numbers are for N = 30. 

71 A n Cn 71 A n Cn 

1 0.8768551 0.6355414 16 14.355486 0.0411725 

2 1.8751352 0.2334135 17 16.420935 0.0468693 

3 2.8725246 0.1380313 18 19.151031 0.0530710 

4 3.8651025 0.0965691 19 22.823289 0.0600117 

5 4.8480195 0.0734911 20 27.880159 0.0679791 

6 5.8147181 0.0585680 21 35.055186 0.0773592 

7 6.7563510 0.0479884 22 45.625953 0.0887028 

8 7.6611613 0.0397497 23 61.958086 0.1028391 

9 8.5126922 0.0328696 24 88.782297 0.1210893 

10 9.2858252 0.0264066 25 136.54610 0.1457075 

11 9.9596476 0.0218315 26 231.60502 0.1808932 

12 10.652940 0.0246335 27 454.24227 0.2354904 

13 11.565606 0.0303418 28 1131.9155 0.3319354 

14 12.125054 0.0000002 29 4536.2877 0.5494628 

15 12.775746 0.0357576 30 81104.991 1.5494668 



FIG. 1. Temperature variation of the dispersion and ab- 
sorption of the linear magnetic susceptibility measured at 
fixed frequencies in fields applied parallel to the c-axis of an 
YBCO-film : (a) B = 0.4 T and (b) B = 4T. 

FIG. 2. Eigcnfunctions f n (r) of the eigenvalue problem, 
Eq. (4), determining the normalized current density J(r, u>) of 
a circular disk. 

FIG. 3. Frequency dependence (a) of the phase angle P and 
(b) of the modulus of the dynamic conductivity cr(u)) evalu- 
ated from the dynamic susceptibility x( w ) measured at 0.4 T. 
Full lines are drawn as guides to the eye. Note that the limit 
of thermally activated flux flow, P = and a(uj) = const., is 
reached only very close to T c . 
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FIG. 4. Dynamical scaling of the linear conductivity, (a) 
Phase angle and (b) modulus of a, both showing the crossover 
from pure ohmic to full screening behaviour at two different 
magnetic fields. 



FIG. 5. Field-effect on the dynamical scaling functions for 
(a) the phase angle and (b) the modulus of (t(uj). The collapse 
of data is achieved by introducing appropiate powers of B in 
the scaling variables. The inset shows the glass-line (Eq. (7)), 
relative to the mean-field transition B C 2- 
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